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I discuss the current status of the comparison between quarkonium-production measurements and 
the predictions of the NRQCD factorization formalism. 
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I. NRQCD FACTORIZATION 

In calculating rates for quarkonium production or de- 
cay, one would like to separate the short-distance heavy- 
quark-antiquark (QQ) annihilation or production pro- 
cess, which has a typical momentum scale p ^ mg and 
can be treated perturbatively, from the long-distance 
quarkonium dynamics, which have a typical momentum 
scale p ^ rnqv and are nonperturbative in character. 
Here mg is the heavy-quark mass, and v is the relative ve- 
locity of the Q and Q in the quarkonium rest frame. Such 
a separation of long- and short-distance scales in heavy- 
quarkonium production and decay can be expressed ele- 
gantly in terms of the effective field theory Nonrelativistic 
QCD (NRQCD) [1]. 

In the cases of inclusive quarkonium production at 
large transverse momentum px and at large CM-frame 
momentum p* , the cross section can be written in a fac- 
torized form as a sum of products of NRQCD matrix 
elements and short-distance coefficients: 



(1) 



The F„(A) are short-distance coefficients. They are, es- 
sentially, the partonic cross sections to make a QQ pair 
convolved with parton distributions. The Fn{A) can be 
calculated as expansions in the strong-coupling constant 
as- The (vacuum) matrix elements involve four-fermion 
operators, which have the form 
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X^k^^(J2\H + X){H + X\)^^k'„X- (2) 
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Here, ip is the Pauli spinor field that annihilates a heavy 
quark, x is the Pauli spinor field that creates a heavy 
antiquark, and k contains Pauli matrices, color matri- 
ces, and covariant derivatives. The operator creates a 
QQ pair in a state with certain color, spin, and orbital- 
angular-momentum quantum numbers, projects it onto 
an intermediate state containing a heavy quarkonium H 
plus anything, and annihilates a QQ pair with specific 
quantum numbers from that state. The operator matrix 



elements contain all of the long-distance (nonperturba- 
tive) physics. They are the probabilities for a QQ pair 
to evolve into a heavy quarkonium. 

A similar factorization formula applies to inclusive 
quarkonium decays. The decay matrix elements are the 
crossed versions of quarkonium production matrix ele- 
ments. Only the color-singlet production and decay ma- 
trix elements are simply related. 

The NRQCD factorization formalism gains much of 
its predictive power from the fact that the nonpertur- 
bative operator matrix elements are universal, i.e., pro- 
cess independent. Although some decay matrix elements 
have been computed on the lattice [2, 3], in general, the 
matrix elements must be extracted phenomenologically. 
The consistency of the phenomenological matrix elements 
from process to process is a key test of the NRQCD fac- 
torization formalism. 

NRQCD also predicts velocity-scaling rules [1], which 
give the leading power behavior of the matrix elements 
as functions of w. It follows that the sum over operator 
matrix elements in Eq. (1) is actually an expansion in 
powers of u. For charmonium, « 0.3; for bottomonium 
?;2«0.1. 

An important feature of the NRQCD factorization for- 
malism is that quarkonium decay and production occur 
through color-octet, as well as color-singlet, QQ states. 
If one drops all of the color-octet contributions, then 
the result is the color-singlet model (CSM). In contrast, 
NRQCD factorization is not a model. It sometimes is 
called, erroneously, "the color-octet model," but it is, 
rather, a consequence of QCD in the limit m, px 3> 
Aqcd- 

A proof of NRQCD factorization would rely both on 
NRQCD itself and on the machinery that is used to prove 
factorization of hard-scattering processes in QCD. Al- 
though it is widely believed that the NRQCD factoriza- 
tion formula (1) can be established by standard methods, 
no detailed proof exists in the literature. It is expected 
that corrections to the factorization formula are of order 
Aqcd/Pt ^'^^ unpolarized cross sections and Aqcd/pt 
for polarized cross sections [4]. 



II. SOME SUCCESSES OF THE NRQCD 
FACTORIZATION FORMALISM 
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Predictions of the NRQCD factorization formalism 
have been confirmed in a number of decay and produc- 
tion processes. In this talk, I will focus on quarkonitmr 
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FIG. 1: The J/i/i cross section as a function of pr- The data 
points are from the CDF measurement [5]. The sohd curve 
is the NRQCD factorization fit to the data given in Ref. [6]. 
The other curves give various contributions to the NRQCD 
factorization fit. From Ref. [61. 



FIG. 2: The cross section for 77 J/ip + X at LEP. The 
sohd curves correspond to the MRST98LO [7] parton distri- 
butions, and the dashed curves correspond to the CTEQ5L [8] 
parton distributions. The upper set of curves is the NRQCD 
factorization prediction, and the lower set of curves is the 
color-singlet-model prediction. The data points are from the 
Delphi measurement [9]. From Ref. [10]. 



production. 



A. Quarkonium production at the tevatron 

The CDF data for J/ ij) production at the Tevatron are 
shown in Fig. 1, along with a fit based on NRQCD factor- 
ization. As can be seen, the color-singlet contributions, 
whose normalizations are known from decay processes, 
undershoot the data by more than an order of magnitude. 
The color-octet contributions offer a possible explanation 
of the data. However, it should be remembered that the 
color-octet matrix elements are obtained by fitting to the 
Tevatron data. A precise test of NRQCD factorization 
occurs only when one uses the fitted values of the matrix 
elements to make predictions for other processes. The 
shape of the data is consistent with NRQCD factoriza- 
tion, although there is a good deal of freedom to change 
the theoretical shape by adjusting the values of the color- 
octet matrix elements. The Tevatron data for ij:' and T 
production are also fit well by the NRQCD factorization 
expressions. 



B. 77 J/V + X at LEP 

As can be seen in Fig. 2, in the case of J /ip production 
in 77 collisions at LEP, comparison of theory [10] with 
Delphi data clearly favors NRQCD over the color-singlet 
model. Theoretical uncertainties arise mainly from un- 
certainties in the color-octet matrix elements and from 
the choices of renormalization and factorization scales. 
The latter uncertainties were estimated by varying the 



scales up and down by a factor of two. The uncertain- 
ties in the color-octet matrix elements are exacerbated 
by the fact that different linear combinations of matrix 
elements appear in J/t/j production at LEP than in J /tp 
production at the Tevatron. 



C. Quarkonium production in DIS at HERA 

The leading-order NRQCD factorization and CSM pre- 
dictions of Ref. [11] for the J/ij] inclusive production cross 
section da/dp^ in ep deep-inelastic scattering (DIS) are 
shown in Fig. 3, along with the HI data. As can be seen, 
the data clearly favor the NRQCD factorization predic- 
tion over the CSM prediction. A similar situation holds 
for da /dQ^. Surprisingly, the cross section da/dz, which 
is differential in the energy fraction (inelasticity) z, is 
not fit well by either the NRQCD factorization or CSM 
predictions. The former overshoots the data at large z 
and undershoots the data small z, while the latter un- 
dershoots the data at all z. It is worth noting that the 
calculation of Ref. [11] disagrees with a number of previ- 
ous results [13-17], which themselves are not fully con- 
sistent. Those discrepancies have not yet been resolved 
completely. 



III. SOME PROBLEMATIC COMPARISONS 
WITH EXPERIMENT 

There are several notable processes for which the com- 
parison of NRQCD factorization predictions with the 
data is less than satisfactory. 
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FIG. 3: J/t/' inclusive production in ep DIS. The vertical axis 
is da/dp^ in units of pb/GeV^; the horizontal axis is the lab- 
oratory in units of GeV'^. The data points are from the HI 
Collaboration [12]. The upper histogram is the leading-order 
NRQCD factorization prediction, and the lower histogram is 
the GSM prediction. Theoretical error bands include esti- 
mates of uncertainties from renormalization and factorization 
scales, parton distributions, and color-octet matrix elements. 
From Ref. [11]. 



A. Polarization of quarkonium at the Tevatron 

The polarization of quarkonium produced at large pt 
at the Tevatron provides a potentially definitive test of 
the color-octet mechanism. Quarkonium production at 
large-pT {pt ^ 4mc for the J/ip) is dominated by gluon 
fragmentation into the quarkonium through the ^Si 
color-octet matrix element. At large p^, the fragmenting 
gluon is nearly on its mass shell, and, so, is nearly trans- 
versely polarized. According to the NRQCD velocity- 
scaling rules, spin-flip interactions are suppressed rela- 
tive to non-spin-flip interactions. Therefore, it is ex- 
pected that most of the gluon's polarization is transferred 
to the J/4> [18]. Radiative corrections and color-singlet 
production dilute this polarization [19-21]. Nevertheless, 
substantial polarization is expected at large px- In the 
J/4> case, feeddown from the tp' and the Xc states is im- 
portant and has now been taken into account [22]. The 
NRQCD factorization prediction for the J/tp polariza- 
tion as a function of px is shown, along with the CDF 
data, in Fig. 4. The quantity a parametrizes the angular 
distribution of the decay leptons in the J/ i/j rest frame: 
da/d{cos6) oc 1 -I- acos^ 9. Here, 9 is the angle between 
the three-momentum of the positive lepton in the J/tp 
rest frame and the boost vector from the J/tp rest frame 
to the CM frame of the colliding hadrons. a = 1 corre- 
sponds to 100% transverse polarization; a = —1 corre- 
sponds to 100% longitudinal polarization. Polarization of 
produced ip' mesons is simpler theoretically, since feed- 
down is not important, but, unfortunately, the experi- 
mental statistics are much poorer. 

As can been seen, the observed J/ip polarization is 
generally smaller than the prediction and seems to trend 
in the wrong direction, decreasing with increasing px- 
However, the experimental error bars are large, and only 
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FIG. 4: J/i/) polarization at the Tevatron. The band is the 
total NRQCD factorization prediction. The other curves give 
the contributions from feeddown from higher charmonium 
states. The data points are from the CDF measurement [23]. 
From Ref. [22]. 



the last data point truly disagrees with the prediction. 
Furthermore, there are large uncertainties in the theoret- 
ical predictions. There are uncertainties in the NRQCD 
matrix elements, which are reflected in the prediction 
band in Fig. 4. There are uncertainties from uncalcu- 
lated contributions of higher order in as, including ef- 
fects from multiple soft-gluon emissions and k-factors. 
Since the polarization depends on a ratio of matrix el- 
ements, it probably is not strongly affected by correc- 
tions to the matrix-elements fits from such higher-order 
effects. Next-to-leading-order corrections have been cal- 
culated for '^Si color-octet fragmentation [19, 24, 25], 
which gives the bulk of the polarization. Corrections to 
the non-fragmentation process could conceivably increase 
the unpolarized contribution by a factor of two. There 
are also large order-w^ corrections to gluon fragmentation 
to quarkonium [26] , but the principal effect of these is to 
change the size of the corresponding color-octet matrix 
element in fits to the Tevatron data, which does not affect 
the polarization prediction. The large order-ti^ correc- 
tions do, however, raise questions as to the convergence 
of the NRQCD v expansion. 

Existing calculations assume that 100% of the QQ po- 
larization is transferred to the quarkonium. Spin-flip cor- 
rections are suppressed only by , not , relative to the 
non-flip part [1]. It could happen that the spin-flip cor- 
rections are anomalously large and depolarize the pro- 
duced quarkonium. It has also been suggested that the 
velocity-scaling rules may need to be modified for the 
charmonium system [27, 28]. These issues should be re- 
solved by a lattice calculation that is in progress [29]. 



B. Inelastic photoproduction at HERA 

Theoretical calculations of the cross section for inelas- 
tic photoproduction of quarkonium at HERA have been 
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FIG. 5: The rate for inelastic quarkonium photo-production 
at HERA as a function of the energy fraction z. The curves 
give total, direct, and resolved color-octet (CO) and color- 
singlet (CS) contributions. The band shows the uncertainty 
in the total contribution that arises from the uncertainties in 
the color-octet matrix elements. The data points are from the 
HI and Zeus results of Refs. [37, 38]. The more recent HI and 
Zeus results of Refs. [39, 40] show a similar behavior. From 
Ref. [6]. 



carried out in the NRQCD factorization formalism by 
several groups [30-36]. The compilation of predictions 
from Ref. [6] and the HI and Zeus data are shown in 
Fig. 5, plotted as function of the energy fraction z. As 
can be seen, the color-octet contribution is poorly de- 
termined, owing to large uncertainties in the color-octet 
matrix elements. Even so, there is little room for a color- 
octet contribution. Furthermore, as is shown in Fig. 6, 
corrections of next-to-leading order in (NLO) [41, 42] 
increase the color-singlet piece by about a factor of two 
at large z and are, by themselves, in good agreement with 
the data. The data shown in Figs. 5 and 6 are for the cut 
Pt > 1 GeV. One can question whether factorization is 
valid at such small values oipr- However, the data differ- 
ential in Pt are compatible with NLO color-singlet pro- 
duction alone at large [6]. It should be noted, though, 
that there are large uncertainties in the NLO color-singlet 
contribution, which arise primarily from uncertainties in 
•nic and cks- The true color-singlet contribution could be 
lower than the central value by about a factor of two, 
leaving more room for a color-octet contribution. 

Near z = 1, the leading-order color-octet contribu- 
tion grows rapidly, in apparent disagreement with the 
data. However, in this, region soft-gluon emission leads 
to large logarithms of 1 — z and also to large corrections 
of higher order in v, both of which must be resummed. 
The resummation of the corrections of higher order in v 
leads to a nonperturbative "shape function" [43]. Both 
the shape function and the resummed logarithmic correc- 
tions significantly smear out the color-octet contribution 



FIG. 6: The rate for inelastic quarkonium photo-production 
at HERA as a function of z, compared with the NLO color- 
singlet model. The band shows the uncertainties that arise 
from as and m. The data points are from the HI and Zeus 
results of Refs. [37, 38]. From Ref. [6]. 
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FIG. 7: HI photoproduction data [39] for da/dz compared 
with the leading-order color-octet (CO), color-singlet (CS), 
and total (CO+CS) contributions for various values of the 
parameter A in the color-octet nonperturbative shape func- 
tion. From Ref. [39]. 



near z = \. The effects of a model calculation of the 
shape function [44] on the leading-order NRQCD factor- 
ization predictions are shown in Fig. 7, along with the HI 
data. The Zeus data [40] show a similar behavior. As can 
be seen, the inclusion of shape-function effects may lead 
to a considerable improvement in the agreement of the 
NRQCD factorization predictions with the data. The re- 
summation of logarithms is expected to have a compara- 
ble effect. (Note that the higher-p^ data shown in Fig. 7 
are more compatible with a color-octet contribution than 
the data shown in Fig. 5). Effects from resummation 
of logarithms of 1 — z and model shape functions have 
been calculated for the process e+e^ J /tp + X [45]. 
For shape functions that satisfy the velocity-scaling rules, 
these effects are comparable in size. It may be possible 
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to use this rcsummcd theoretical prediction to extract 
the shape function from the Belle and BaBar data for 
e+e" J/ip + X and then use it to make firm predic- 
tions for J/'0 photoproduction near z = 1. 

C. Double cc production at Belle 

For the exclusive double charmonium process e+e^ 
J/ip+T]c, the Belle Collaboration measures a cross section 
times a branching ratio into at least four charged tracks 
of 46 ±6^9 fb [46]. In contrast, leading-order calculations 
predict a cross section of 2.31 ± 1.09 fb [47-49]. There 
are some uncertainties from uncalculated corrections of 
higher-order in and v and from NRQCD matrix el- 
ements. However, because this is an exclusive process, 
only color-singlet matrix elements enter, and these are 
fairly well determined from the decays J/^ — > e+e~ and 
Vc 77- 

Since the Belle mass resolution is 110 MeV but the 
J/i'-Vc mass difference is only 120 MeV, it has been sug- 
gested that some of the J/ip + rjc data sample may consist 
of J/tp + J/tp events [50, 51]. The state J ftp + J/tp has 
charge-parity C = +1, and consequently, is produced in 
a two-photon process, whose rate is suppressed by a fac- 
tor (a/tts)^ relative to the rate for J/ip + rjc- However, as 
was pointed out in Refs. [50, 51], the two-photon process 
contains photon-fragmentation contributions that are en- 
hanced by factors (-Ebcam/2TOc)'^ from photon propaga- 
tors and log[8(i?bcam/2mc)^] from a would-be coUinear 
divergence. As a result, the predicted cross section 
cr(e+e- J/ij) + J/ip) = 8.70 ± 2.94 fb is larger than 
the predicted cross section (j{e^e~ — » J/ip + rjc) = 
2.31 ± 1.09 fb, although corrections of higher order in 
a and v likely reduce the former prediction by about a 
factor of three. These predictions spurred a re-analysis of 
the Belle data [52] , with the result that there is no signifi- 
cant J/tp + J/tp signal observed [(T(e+e^ J/ip + J/ip) < 
7fb]. 

There are also Belle results on inclusive double- 
charmonium production. For the ratio Rj/^ = 
a{e'^e~ ^ J/tp + cc)/a{e^e~ J/tp + X), the most re- 
cent Belle analysis yields Rj/^ = 0.82 ± 0.15 ± 0.14, with 
Rj/^ > 0.48 (90% confidence level) [46]. Predictions 
based on NRQCD factorization [53 55] give Rj/^^ ~ 0.1. 
The measured and predicted J/tp + cc absolute cross sec- 
tions also disagree by almost an order of magnitude, with 
the Belle result [56] being about 0.6 1.1 pb and the pre- 
diction [53-55] being about 0.10-0.15 pb. This prediction 
is based only on the color-singlet contribution. However, 
corrections of higher order in v, including color-octet con- 
tributions are not expected to be large. Neither are cor- 
rections of higher order in as ■ 

The discrepancies in the double cc inclusive and exclu- 
sive cross sections are among the largest in the standard 
model. Theory and experiment differ by about an order 
of magnitude — a discrepancy which is larger than any 
known QCD k-factor. It is important to recognize that 



these discrepancies arc problems not just for NRQCD fac- 
torization, but for perturbative QCD (pQCD) in general. 
In the case of the cross section for e+e^ J/tp + rjc, one 
obtains exactly the same predictions in the NRQCD fac- 
torization [47, 48] and light-front-QCD [49] formalisms. 
With regard to the fraction of J/tp + X events that are 
J/tp + cc, it is difficult to see how any perturbative cal- 
culation could give a value as large as 80%. The color- 
evaporation model, for example, proceeds through the 
same Feynman diagrams as the NRQCD factorization 
calculation, differing only in the specific treatment of 
the evolution of the cc pair into quarkonium, and would, 
therefore, be expected to give a prediction that is not too 
different from that of NRQCD factorization. 

Clearly, it is very important to have independent 
checks of the Belle double-charmonium results, such as 
could be provided by the BaBar collaboration. If the 
Belle results are confirmed, then we would be forced to 
entertain some unorthodox possibilities: there are new 
charmonium production mechanisms within the standard 
model that have not yet been recognized, pQCD is in- 
applicable to double-charmonium production, or physics 
beyond the standard model plays an important role. It 
would be very surprising, however, if the last possibility 
could manifest itself at such low energies. 

IV. SUMMARY 

The NRQCD factorization approach provides a sys- 
tematic method for calculating quarkonium decay and 
production rates a double expansion in powers of as 
and V. Calculation of production rates also relies upon 
hard-scattering factorization, for which corrections are 
suppressed by powers of Aqcd/pt- NRQCD factoriza- 
tion has enjoyed a number of successes, for example, in 
quarkonium production at the Tevatron, 77 ^ J/tp + X 
at LEP, and quarkonium production in DIS at HERA. 
Other experimental tests are, so far, more problematic. 
These include qiiarkonium polarization at the Tevatron, 
inelastic quarkonium photoproduction at HERA, and 
double cc production at Belle. The Bclk; double cc pro- 
duction results present a severe challenge to pQCD. It 
would be very useful for the BaBar collaboration to check 
these results. In other cases, inclusion of corrections 
of higher order in and v and resummation of soft- 
gluon effects and endpoint corrections of higher order in v 
should help to achieve agreement between theory and ex- 
periment. More precise theoretical predictions are ham- 
pered by uncertainties in the NRQCD matrix elements. 
Lattice calculations can help to pin down the decay ma- 
trix elements, but it is not yet known how to formu- 
late the calculation of production matrix elements on the 
lattice. This is an exciting time for heavy-quarkonium 
physics, with a great deal of experimental and theoret- 
ical activity in quarkonium decay and spectroscopy, as 
well as production, and we can expect to see continuing 
progress on the many challenging problems that remain. 
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